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Abstract 

q\ We analyze isospin breaking through quark mass differences and virtual photons in the pion-nucleon scattering lengths in all 
C^) physical channels in the framework of covariant baryon chiral perturbation theory. 
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Isospin violation in the Standard Model is driven by strong 
and electromagnetic interactions, that is by the differences in 
the light quark masses and charges, respectively. As already 
stressed by Weinberg, the pion-nucleon scattering lengths of- 
fer a particularly good testing ground for strong isospin vio- 
lation |1]. This problem was addressed in the framework of 
heavy-baryon chiral perturbation theory (ChPT) in a series of 
papers about a decade ago JHHHillsll^t]. Recently, new interest 
arose in high-precision calculations of the pion-nucleon scatter- 
ing lengths. First, the accurate measurements of the character- 
istics of pionic hydrogen and deuterium allow one in principle 
to extract certain jiN scattering lengths to high precision. This, 
however, is only possible if isospin breaking is taken into ac- 
count consistently. In the case of the strong energy shift of the 
ground state of pionic hydrogen one needs the isospin-violating 
contributions to a„- p ^„- p . In [7], these have been determined 
at third order in the chiral expansion, 0(p 3 ), in a covariantly 
regularized form of baryon ChPT |8||. [For a recent review on 
baryon ChPT, see [9|].] As for the width of pionic hydrogen, 
the knowledge of the isospin-breaking corrections to a n - p ^ n 
is required. In the analysis of pionic deuterium isospin viola- 
tion is particularly important, since the nd scattering length at 
leading order is proportional to the small isoscalar scattering 
length a + and therefore chirally suppressed (cf. lUOIO . Since 
Re Qjrd cc a n -p^- p + a^n^n-n + few-body corrections, we may 
improve at least the two-body contributions by extending the 
isospin-breaking corrections to a„-n-Hrn to 0(p i ). Second, as 
has been stressed in particular by Bernstein, threshold pion pho- 
toproduction offers the unique possibility of measuring the so 
far undetermined n°p scattering length and gives access to the 
charge exchange scattering length a n +„^ n o p , see 111 111 and the re- 
cent review 11211 . Such measurements are becoming feasible at 
HlyS and at MAMI. In view of these developments, it is timely 
to extend the work of [7] to all charge channels in pion-nucleon 
scattering. 



2. Formalism 

We start the description of various formal aspects of jiN scat- 
tering at threshold with the kinematics. The momenta of the 
nucleon and pion in the initial (final) state will be denoted by 
p (p') and q (q'), respectively, their masses by nii (ntf) and M\ 
(Mi). m p , m n , M„, and M n o, are the masses of proton, neutron, 
and charged and neutral pion. We define the isospin limit by 
the charged particle masses m p and M„. Working at first order 
in isospin breaking, i.e. at 0(e 2 , m& - m u ) = 0(6), we only need 
contributions linear in A n = M% - M^ and An = m a - m p . 

For elastic scattering, the kinematics at threshold are deter- 
mined by 



s = (nti + Mi) 



— q=— q', t = 0. (1) 

Mi Mi 



([T]) is modified for the charge exchange reactions (cex) accord- 
ing to 

p*p\ q*q', t = -An + Mn (m] - m\ + h n ). (2) 

nii + M„ 

In loop contributions that only start at 0(p 3 ), these kinematical 
relations may be chirally expanded, leading to 

■ffi 

s = 0« p + M„f, t = -A*, p=p' = -j-q. (3) 

Note that still q' must not be replaced by q, since the difference 
is of the same chiral order as q and q' themselves. 

The pion-nucleon scattering amplitude T„n is parameterized 
in terms of the two amplitudes D(s, f) and B(s, f) according to 

T„ N = u( P ')(D(s,t) - -[q',q}B(s,t)\u(p), 

\ 2(n%i + ntf) I 



nif) 

u(p')u(p') - 2mi, u(p)u(p) = 2m;. 
In the isospin limit, T„n may be decomposed as 



nab 



T + 6 ab + T--[T a ,T b ], 
2 



(4) 



(5) 
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where a (b) is the isospin index of the outgoing (incoming) pion 
and t' are the Pauli-matrices. Using the Condon-Shortley phase 
convention, the physical amplitudes are related to T + and T by 



T = 


T 


— T-x+n — T n+n^>n + n 


= T + + T , 




T + = 

* 7I + p — 


Tjr + p— >7T + p 


- T = T 


= T + - T~, 




ycex _ 

1 7T-p - 




ycex _ y 
- 1 n + n - 1 n + n-*7t°p 


= -yf2T~, 




Tjflp = 


Tjfi> p^>7f> p 




= T + . 


(6) 



For the elastic channels only D(s, t) contributes at threshold, 
whereas we find for the charge exchange reactions 
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2 y/m n m p 



1 + 
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Am 2 p 



B 
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(7) 



where and B± r denote the amplitudes evaluated at thresh- 
old. The correction factor in front of D stems from the expan- 
sion of the Dirac spinors around the isospin limit. Since the 
prefactor is already of first order in 6, B^ may be evaluated 
assuming isospin symmetry to relate it to isovector threshold 
parameters 11311 . 



fl 0+ 

B thr = 87rm p[ ^2 
1 



1 + 



= 2^2(1 +4m p c 4 ) + 0(p), 



B thr = - V2B thr , (8) 



where a, ± denotes the isovector scattering lengths with orbital 
momentum I and total angular momentum / ± i. For brevity, 
we will use a* = a? for the S-wave isoscalar and isovector 
scattering lengths. Equation ([8]l also shows the leading chiral 
representation of B~L. All relevant terms of the effective chiral 
Lagrangians defining the corresponding low-energy constants 
are collected in Appendix IA1 

The S-wave scattering length a for elastic scattering of scalar 
particles is related to the amplitude T(s, t) by 



a - 



8ttVs 



T{s,t) 



IpHO 



(9) 



where |p| is the center-of-mass momentum. This result is gen- 
eralized to pion-nucleon scattering by 



^elastic — 



-D 



.elastic 



4^(mi + Mi) thr ' 



4n(rm + M n ) 



(10) 



The isospin-symmetric contributions to the scattering lengths 
have already been worked out in |fl4n . Adapted to our notation 
they read 



m p M 2 n 



An(m p + M n )Fl 

rripMx 
87T(m p + M n )Fl 



I _ A. + 2(c 2 + c 3 - 2ci) + 
I Am p 

( 1 + i^ + ^ (l - log 



1 - log -f 

V 2 



+ 8M 2 n (d\ + 4+4 + 24) + 



2M. 



(11) 



As soon as we take into account virtual photons, we have to 
specify more carefully what we mean by D± r due to the ap- 
pearance of threshold divergences. First of all, we subtract all 
one-photon-reducible diagrams, since they diverge ~ l/t, and 
denote the result by D. The additional divergences due to pho- 
ton loops may be regularized in the form 

UQ^mQ( s , t ))\ =^ +y62 log^ + Ahr + tf(|p|), (12) 
v ;| IpH0 Ipl fi c 

where a = e 2 /An denotes the fine structure constant, #c(lpl) the 
infrared divergent Coulomb phase given by 



Uc m v 

0c(|p|) = -—log — , 
V Ipl g 2|p|' 



(13) 



and p c = m p M n /(m p + M n ) the reduced mass of the incoming 
particles. Q accounts for the charges of the particles involved, 
explicitly 



&rp = -2, 0^=2, Q™ x p = -1, 



(14) 



and Q - for the remaining channels. For consistency rea- 
sons, the contribution from B(s, t) to the charge exchange re- 
action should be multiplied by the same phase as D(s, t). Since 
af9c(|p|) is of first order in isospin breaking, this does not lead to 
additional terms at the order considered here, but makes it ob- 
vious that r?c(|p|) drops out of physical observables. The coeffi- 
cients f3i may be related to the scattering lengths a. The explicit 
relation 

Px=- n -Qe 2 M K a (15) 

is confirmed by our calculation at leading order in the chiral 
expansion, but can be proven to hold in general in the frame- 
work of non-relativistic effective field theories [15. 

d. The 

coefficient /?2 only appears at two-loop level. 

3. Analytic results 

The topologies of the Feynman diagrams contributing at 
threshold are displayed in Fig. [TJ There are significantly less 
diagrams than above threshold, since many diagrams which are 
formally of 0(p y ) start only at 0(p ) for the following reasons. 
Firstly, the leading term of a particular diagram can be deter- 
mined by simplifying the numerators according to chiral power 
counting, making use of (HJ and J3). With X = p + q and loop 
momentum k a typical example for such a simplification is 

M 

(2 - t + m p )qy 5 {f -1 + m p ) -> (jj> + m p )—py 5 (p + m p ) 

Win 



= -75 H>(-# + m p ){p + m p ) = 0. 



(16) 



Secondly, all s-channel one-particle-reducible diagrams involve 
structures of the type (2 + m p )fysu(p), whose leading part van- 
ishes at threshold, since 



M„ 

(t + m p )qysu{p) -> {f + m p )—fysu(p) 

trip 

= —pli(-f + m p )u(p) = 0. 
m p 



(17) 
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Figure 1: Loop diagrams for nN scattering at threshold. Solid, dashed, and 
wiggly lines, denote nucleons, pions, and photons, respectively. Crossed dia- 
grams and diagrams contributing via wave function renormalization only are 
not shown. 



The M-channel diagrams are treated analogously. Unfortunately, 
both arguments only work for q and not for q' in the charge 
exchange reactions, unless the diagram in question is already 
of order 0(5); but eventually one can show that all diagrams 
which may be omitted in the case of the elastic channels do not 
contribute to the charge exchange reactions either. 

Concentrating on the analysis of the isospin-breaking shifts 
in the scattering lengths, we obtain the following results for the 
reaction channels on the proton (the neutron channels can be 
found in AppendixlBli: 



Aa K - p = 
Art ' = 



*n p 



(a + + a ) = Aa + + Aa + i Im a„- 



(a + — a ) = Aa + 



Aa 



Aa - - 



m p \4A n e 2 
4n(mp + M n ) { F* 2 

8^(33^, 2 \ 
32jiF 2 \4F 2 . e ), 
m p M„ 



4n(m p + M„) { 32n 2 F A K 



Ml 

3+log-f 



+ 8A, ^fL 
+ Fi d5+ IG^n 



Ml 

1 + 41og2 + 31og-f 

y" 2 



\ma K 



Mi 



, - , — ^- VA^ - 2M n A N + £ 8aM " 

" 4n(m p + M n )\%nF* V N 4nF% 



^ u jrp u irp T 



V5< 



V2; 



4n(m p + M„) { 2 



e 2 fi 



g\A n M„A 



4F 2 n m p 



8MA 

F 2 d 5 



\92n 2 Fi 

e 2 M., 



M„A X ( i \ 

+ ^f\ B ^-4-Fl) + 

I M 2 \ 
-^-7$*+p-5«*)log-Sj 



Ml 

3*m \ 5 + 3]og ^ 



M„A^ 
4F 2 n m p 



(i + 2«i) 



e M„ t , r . r 20, , rs \] 
^\F 2 ^ 1 -2k\ + k\ + -{k\ + k\)^, 



Aa n o 



A, 



M 2 " + 47T(m p + M n ) \ Yl%nF* 



3g\M n A n 



Ml y/A„ + 2M K A N 2c 5 B(m d - m u ) 



Fl 



(18) 



We wish to point explicitly to the square-root-like terms in 
Ima„- P and Aa,fl p , which are caused by the unitarity cusps 
due to the different thresholds of the n°n and n + n intermedi- 
ate states, respectively. These cusps can be calculated exactly 
at threshold, which we will illustrate for the imaginary part in 
Sect. 14.31 Since the cusp is of order 0( yfd) and thus enhanced 
compared to the other isospin-breaking effects, we also take 
into account the correction by An, although this is formally an 
0(p 4 ) effect. Nevertheless, it contributes ~ 30 % to the differ- 
ence between a n o p and a n o„ (see AppendixlBt. 

We have performed the following checks on our calculation: 
the amplitudes are ultraviolet-finite, all ultraviolet divergences 
due to loops are canceled by the infinite parts of the countert- 
erms (as calculated in |7[]). Thus, only the renormalized coun- 
terterms appear in ( fT8l . They compensate the scale dependence 
generated by the chiral logarithms, such that the final results are 
independent of the renormalization scale fi. Furthermore, the 
infrared divergences caused by virtual photons cancel among 
themselves, as they should. 

A useful way to quantify isospin-breaking corrections in 
terms of measurable quantities is the so-called triangle relation 
that vanishes in the isospin limit. It is defined as 



R = 2 



2 R + p - a„- p — V2a" 



cex 
P 



"bfip - a„- p + V2a^ 



(19) 



where only the real parts of the scattering lengths are inserted. 
At first order in 6 we obtain 



R = 



e 2 fi g\^n M n A 



(1+2 i 

4n{m p + M„)a-\ 2 ' 4F 2 m p 4ffm^ + 8a) 

3M n A n M„A„ M K A„ I M 2 

+ — — —B thr - - ^ 4 + log — f 



\6F 2 n m 2 p 



4m 



2 thi 



4%n 2 Fi 



g\M n A n l M 2 n \ e 2 M n I M 2 A 

e 2 g\M n 



16 ^ M l + 41og2 + 31og^ 

p 2 M p 2 M 



(20) 



We refrain from constructing an isoscalar triangle relation from 

I — I 

the three elastic pion-proton scattering lengths (cf. R\ in yfl); 
such a relation can easily be read off from the results in ( fT8l . 
It depends on the low-energy constants f\, fa, and c\, and, as 
we will see in the following section, therefore cannot be very 
well constrained, such that no additional information beyond 
the shifts in the individual scattering lengths is provided. 



4. Numerical results 



4.1. Low-energy constants 



The most precise values for a + and a stem from an analysis 
of pionic hydrogen and pionic deuterium data lioll 



4.2. Triangle relation, scattering lengths 

The triangle relation R can be determined rather well since 
fx, the 0{p 2 ) LEC which is least known, drops out. The cen- 
tral value is obtained by inserting the above LECs and (ETT l into 
(f2Qb . As for the error, we are faced with the following combina- 
tion of electromagnetic LECs whose uncertainty is not known: 



a + = (1.5±2.2)-10 _J M; 1 , a = (85.2 ±1.8)- KTM^ 1 . (21) 

In addition, the authors extract the electromagnetic low-energy 
constant (LEC) f\ = -2.1^1 GeV -1 . f and C5 can be deduced 
from the mass difference between proton and neutron. This 
mass difference comprises electromagnetic as well as strong 
contributions 



m n -m p - -4Bc 5 (m d - m u ) + fe 2 F 2 n , 



(22) 



which may be disentangled by means of the Cottingham for- 
mula 1 17]. The result of this procedure is f = -(0.97 ± 
0.38) GeV Bc 5 (m d - m a ) = -(0.51 ± 0.08) MeV. 

In II 1 811 - various previous analyses of c\ are briefly reviewed 



and combined to c\ = -O^+^GeV 
F*dZ(M) = (0.6 ± 3.0) • 10~ 3 



For d r 5 , we will use 



s , f , , - _ , ,, , specifying the renormalization 

scale to fi = 1 GeV. The central value is the mean of the values 
quoted in [ 1 3{| (translated to our conventions regarding '), 
where a low-energy theorem linking d' 5 to a certain subthresh- 
old parameter of nN scattering is derived. In the spirit of the 
treatment of c\, we estimate the error by investigating the ef- 
fects of higher orders in this low-energy theorem. Neglecting 
the fourth order contribution would shift F^dl by 



M K 64m pCl + g 2 A [2(4 + g A ) + V2 log(l + V2)] 



16 



32^ m„ 



-3 ■ 10" 



(23) 

(for the central value of ci). The resulting uncertainty ensures 
consistenc y w ith most values for available in the literature |@, 



We now turn to the determination of Z? thr . Values for a x _ and 



(-81 ±6)- 10~ 3 M; 3 (24) 



a[ + can be found in ||2Jl |2J, |23J, |24J] . Using 



-3i*-3 



(-12±2)-10 _J M; 



'1+ 



yields 



B thr = (0.60: 



:0.06)- 10~ 3 MeV~ 2 . 



(25) 



Since the main source for the determination of C4 are nN thresh- 
old data, it seems more reliable to apply the threshold parame- 
ters directly. 

Estimates of the meson-sector electromagnetic LECs are 
given in [25] using resonance saturation [26, 27]. Unfortu- 
nately, this method does not provide reliable error estimates. 
The central values for the in question are k\ = 10.9 • 10 -3 , 
k\ =0.7-10- 3 ,/t'3 = 3.9-10- 3 ,^ = -1.3-10- 3 (all at// = 1 GeV). 
Since g r 6 , g^, and g^ are not known, they are set to zero in the 
numerical work. Particle masses and decay constants are taken 
from il, in particular = 92.2 MeV and |g A | = 1.2695. 



k\-2k\ + F 2 J4gl + g l 7+ 4gl). 



(26) 



Naively one would assign the order-of-magnitude errors 1 /167r 2 
to each LEC and add the individual contributions in quadra- 
ture. However, this may underestimate the uncertainty in case 
the variation of the renormalization scale fi by a factor of 
e = 2.718..., controlled by the corresponding yS-functions, 
induces a change significantly larger than 1/I6n 2 . Assuming 
1 GeV to be a "natural" scale for hadronic processes, this run- 
ning covers the energy range the physics we consider should be 
sensitive to. Estimating the uncertainty by varying the LECs 
according to their /^-functions in a correlated way has the fur- 
ther advantage of being independent under redefinition of the 
Lagrangian. The result of this procedure is 



R = (1.5 + 0.2y 2 ± 0.03 fl - ± 0.03^ + 1.1 LEC ) % 
= (1.5 + 1.1)%, 



(27) 



where the different contributions to the error are denoted by 
a subscript, "LEC" standing here and in the following for the 
corresponding combination of LECs with unknown error. The 
final uncertainty is obtained by adding the individual contri- 
butions in quadrature. Naive order-of-magnitude arguments 
would reduce the error significantly to 0.4 %. The large error 
in (|27| > is dominated by the g x { , as may be seen from their B- 
functions |7| 777 = -9/2 - 2Z(5g 2 A + l)/3 = -9.4, m = -2i] 6 = 
3(4g 2 A - l)/2 + 2Z(5g A + l)/3 = 13.1, which are by no means 
of order 0(1). 

We now turn to the isospin-violating contributions to the in- 
dividual scattering lengths. The procedure as described above 
yields 

Afl ff - p =(-3.4^ : 2 fi ! 3 ^ /i ±0.2 /2 ±0.6, /5 ±1.2 LE c + 5.0/)l0- 3 M; 1 

= (-3.4^; 3 +5.0/)- io^m; 1 , 

tu* p =(- 5.3!^!;^ ± 0.2 /2 ± 0.6 ds ± 1.2 LEC ) • lO^M; 1 

= -5.3!«.10-X 1 . 
Aa™ = (0.4 ± 0.2 /2 ± 0.8 d5 ± 0.04^ ± 0.4 LEC ) ■ lO^M; 1 

= (0.4±0.9)-10" 3 Af; 1 , 
Aaj> p = (- 5.2 ± 0.1 fl+ ± 0.2 C5 ) • lO^M; 1 



(28) 

corresponds to relative 
1+4.9 < 



= (-5.2±0.2)-10" 3 M; 1 . 
Discarding the imaginary part, 

changes compared to the isospin limit of — 3.9*ig% in a„- p , 
+6.4+™ % in a„* p , and (-0.4 ± 0.8) % in a c n ex p . Due to the poor 
knowledge of a + , the corresponding normalization of Aa„o p is 
not very meaningful; note that the isospin-breaking shift Aa„o p 
in d28l is significantly larger than a + . As already pointed out 



4 



in J3l, the large isospin-breaking corrections to the charged- 
pion elastic channels can be traced back to the triangle graph 
(s$) (see Fig.lTJ, which however only yields a rather small con- 
tribution to the charge exchange reaction. In contrast, isospin 
violation in the neutral-pion elastic channel is predominantly 
due to the cusp effect enhanced by V<5. The large uncertain- 
ties in Aa„± p are dominated by f\ and c\ that are part of Aa + in 
( [T8| >. therefore appear in the same combination with a + in both 
channels. 



cex 
P 



4.3. Imaginary parts 

Exact expressions for the imaginary parts of a^-p and a^ e) 
generated by the n°n and yn intermediate states can be obtained 
using Cutkosky rules, expressing the vertices at threshold by 
scattering lengths and electric dipole amplitudes £0+ encoun- 
tered in the context of pion photoproduction. Retaining all chi- 
ral orders, the resulting imaginary parts up to 0(6) are 



1m 



l 7l p 

.cex 



y/mp + M n 
M n E a+ (n-p) 
Op + M n ) 



(29) 



The experimental value for Eq + (tt~p) taken from ||29|l and the 
leading term of its chiral expansion calculated in 113011 up to 
0(ep 3 ) are 



Eo+(n p) 



V2, 



■+0(ep) 



(-31.5+0.8)- io^m; 1 , 



(30) 



whereas E(, + {n°n) only starts at 0(ep 2 ) (explicit expressions are 
given in 113 111 ): therefore both contributions to Ima^ e * in d29l 
are suppressed by at least one chiral order. Unfortunately, 
Eo+(7T°n) is not directly accessible in experiment. Combining 
deuterium data with ChPT predictions [32] yields 

E Q+ (n°n) = (2.1 + 0.5) • lO^M; 1 . (31) 

Inserting the chiral expansions of Eo + (jr~p) and a~ into < f29b 
reproduces the imaginary part of a„- p appearing in 11 1 St , One 
can easily check that the difference between (1291 and its chiral 
expansion is mainly due to An, which justifies our treatment of 
the cusp effect in Sect. [3] Separating strong (first number) and 
electromagnetic contributions, we obtain numerically 



lma„- p = ((2.91 + 0.12) + (1.86 ± 0.09)) ■ lO^M^ 1 



= (4.77 ±0.15)- 10 _3 M;\ 



(32) 



Im<4 e * = ((-0.04 ± 0.05) + (-0.12 ± 0.03)) ■ 10 
= (-0.16 ± 0.06) ■ io^m; 1 . 



3 M-J 



Finally, the above results may be checked based on the ob- 
servation that the ratio 



Imfl,-,, 

* P I strong 

Imfln-nL,, 
" P\EM 



= 1.57 + 0.10 



(33) 



should correspond to the so-called Panofsky ratio P = 
<r(n~p — > ji n)lcr{n~p — > yn). Indeed, its experimental value is 
found to be P = 1.546 ± 0.009 (cf. SB). 



5. Comparison to earlier work 

Our result for n~p — > n~ p agrees with J3]. In J3l, a simi- 
lar analysis of isospin breaking was performed in heavy -baryon 
ChPT, switching off virtual photons. This corresponds to 

„2 



+ 0, An # 0, fie 2 * 0, f 2 e 2 * 0, m d + m u , e L = 0. 



(34) 



Furthermore, isospin-breaking effects due to the Dirac spinors 
are neglected and is expressed by c\. We have checked 
explicitly for the triangle relation, for the charge exchange re- 
actions, and for the neutral-pion elastic channels, that our re- 
sults coincide in this limit. Numerically, we find R = (0.74 + 
0.21) %, which is compatible with the numerical value /?fms = 
(0.9 ... 1.1) % quoted in [3]. Both values slightly differ, since 
the denominator is not expressed by a' (the additional LECs 
needed are taken from [21]) and since the isospin limit is de- 
fined as the average between charged and neutral particles. 

Virtual photons were taken into account in [6] in order to 
study isospin violation above threshold. Unfortunately, a direct 
comparison is not possible, as no analytic expressions for the 
amplitudes are provided. Even more, also a numerical compar- 
ison is difficult due to a conceptual difference: in |]6j], the elec- 
tromagnetic corrections were used to pin down the LECs from 
experimental data, and thereafter applied to extract the strong 
amplitude. In particular, electromagnetic contributions to the 
particle masses were switched off. Thus, the quoted isospin- 
breaking effect of -0.7 % for the triangle relation in the S-wave 
refers to strong isospin violation only. 

6. Summary and outlook 

In this letter, we have systematically analyzed isospin vio- 
lation in the nN scattering lengths in all channels, including a 
detailed estimate of the theoretical uncertainties. The extension 
of this analysis beyond threshold will be the subject of future 
work, to which we also refer for details of the calculation [35]. 

We find that isospin violation is quite small in n~p — > n°n, at 
the order of one percent at most, whereas the charged-pion elas- 
tic channels display more sizeable effects on the few-percent 
level. In particular, the so-called triangle relation that vanishes 
in the isospin limit is violated by about 1.5% consistent with 
earlier findings in heavy -baryon ChPT and inconsistent with the 
5-7% deviation extracted from the data at lowest pion momenta 
in [36, 37^. In addition, we find a substantial isospin-breaking 
correction to the neutral-pion-proton scattering length. In view 
of these results, further experiments e.g. at HlyS and MAMI 
are urgently called for. 
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A. Effective Lagrangians 

We will use the effective Lagrangian for nucleons, pions, and 
virtual photons, as constructed in [7], whereof we actually need 
the following terms: 

= jfrf> + c sr •) + 1 4f + X + 

i=\ i=l i=0 

F 2 

£» ) + 4 e_) + A = -fidfrfdpU + X "U + U f x ) 
+ ZF 4 (QUQU f ) - if^ - i(fl^A") 2 , 

J^ 2p2) = F^uU^UKkiiQ 2 ) + k 2 (QUQU^) 

+ k 3 ((d>'U f QU){d li U f QU) + (d^UQU^id^UQU^)) 
+ k 4 (d f 'U' t QU)(d tl UQU' t )}, 

4f =*{a9-m+igj*7 5 }¥, 



4f ' = ¥{ci<* + > - -^{u^WW + h.c. 

+ -^{U^) + -CA^iu^, Ky] + C 5 ^ + }*P, 

|L[ Mm ,[D v ,^]]DV_^.[ 



^3 



[^, [Z) v , ^]](Z^D v D' 1 + sym) 



12m3 L ^' 
+ ^d 5 \x-,u fl \DK}v + h.c, 
iF : 



(35) 



-4 = 2^*{ft<(2 + ><(2-«^ +^7/8<e ± « A ,>G T £"'}^ + h.c., 

where (A) denotes the trace of a matrix A, A - A - (A)/2 its 
traceless part, T'((9 + h.c.)T' = TCMf + h.c. for an operator O and 

d M U = 8^11 - iA M [Q, U], x = 2B diag(m u , m d ), U = m 2 , 

= a„A y - 8 V A^ Q=|diag(2,-1), Q = ediag(l,0), 

o P = 3^ + r M , = -(« t (5 /1 - /Qa^h + - /gA^y 1 ), 

,Y± — "'A'"' ± ux^ u, — iyi 1 (dft—iQA^u — u{d ll —iQA l ^)u i y 

Q± = ^(uQu 1 + M 1 Qu), [Dp, u v ] = d^Uy + [r,,, u v ]. (36) 

Y = (p, n) T contains the nucleon fields and the matrix U col- 
lects the pion fields in the usual way. F is the pion decay con- 
stant in the chiral limit and is replaced by its physical value F„ 
by means of 

while the chiral-limit axial charge g may be identified with its 
physical value g\, since axial contributions only start at 0(p 2 ) 
at threshold, m denotes the nucleon mass in the chiral limit. 



Changing the version of J from l38ll to ll39ll results in a 
redefinition of dy This d r 5 is related to our d l 5 by 

F 2 n d\{p) = F 2 d ! 5 (p) + I$Qi). (38) 
Note that in this convention f. disappears in ( fTTT ). 

B. Pion-neutron scattering lengths 

The strong contributions to the remaining channels are deter- 
mined by charge symmetry (the discrete subgroup of the gen- 
eral isospin transformations that only exchanges u «-> d on the 
quark level), such that only the electromagnetic parts have to 
be calculated explicitly: the pion mass difference alone cannot 
contribute to charge-symmetry breaking. [How to simplify a 
calculation of isospin-breaking effects by such considerations 
is explained in more detail in B40I1 .1 The results are 

Aa n+ „ = a^„ - (a + +cT) = (- 4.3^ + 6.0/) ■ lO^M; 1 



Aa K - p + 



47r(m p + M„) 



f 2 -2M !r (2g l 6 +gl) 



M 



Aax-n = a„-„ - (a + - a) = -6.2^-3 . lO^M: 1 



= Aa n + n + 



47r(m p 

Aa™ = a~ + V2V = (2.3 ± 0.9) • lO^M; 1 
V2m D f M n A N 



= Ad 



(l+2gi)-e 2 f 2 , 



Aa 



71 n ^rrn 



- Aa 



r ' 4^(mp + Mj\2F> p 

(-1.8 ±0.2)- lO^M^ 1 
m n ( 4c 5 B(m d - m u ) 



71°/! 



47r(m p + M„) 



Fl 



M- 



+ VA^ + 2M^A N - VA* - 2M^A N ) 



(39) 



fl ™«' which is accessible through the cusp in neutral-pion pho- 
toproduction on the proton, receives only moderate isospin- 
breaking corrections ((-1.9 ± 0.8)%), whose uncertainty is 
rather well-controlled. The correction to the two-body contri- 
bution to Re a n d oc 2(a + + A5 + ) + . . . displays the same depen- 
dence on fi and ci as 2(a + + Aa + ). It is determined by 



Act 



g l k M„mA n 2 \\ 

32nFi\ m )]' 



An(m 

,2 



which, in addition to fi and CK_also includes a sizeable fixed 
shift (cf. the discussion in 14 ll 14211 ). Finally, we point out that 
the remnants of the cusp effect contribute roughly one third to 
the difference a^> p - a^ n = (-3.4 ± 0.4) ■ 10~ 3 M n l and hence 
modify Weinberg's prediction [1] significantly. [This is in ap- 
parent contrast to the finding in Jj] where the complete 0(p 4 ) 
corrections to a n o p - a„o n have been calculated; however, the 
result for the cusp is incorrect.] 
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